Geometrical frustration effects on charge-driven quantum phase transitions 
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The interplay of Coulomb repulsion and geometrical frustration on charge-driven quantum phase 
transitions is explored. The ground state phase diagram of an extended Hubbard model on an 
anisotropic triangular lattice relevant to quarter-filled layered organic materials contains homoge- 
neous metal, 'pinball' and three-fold charge ordered metallic phases. The stability of the 'pinball' 
phase occurring for strong Coulomb repulsions is found to be strongly influenced by geometrical 
frustration. A comparison with a spinless model reproduces the transition from the homogeneous 
metallic phase to a pinball liquid, which indicates that the spin correlations should play a much 
smaller role than the charge correlations in the metallic phase close to the charge ordering transi- 
tion. Spin degeneracy is, however, essential to describe the dependence of the system on geometrical 
frustration. Based on finite temperature Lanczos diagonalization we find that the effective Fermi 
temperature scale, T* , of the homogeneous metal vanishes at the quantum phase transition to the 
ordered metallic phase driven by the Coulomb repulsion. Above this temperature scale 'bad' metal- 
lic behavior is found which is robust against geometrical frustration in general. Quantum critical 
phenomena are not found whenever nesting of the Fermi surface is strong, possibly indicating a first 
order transition instead. 'Reentrant' behavior in the phase diagram is encountered whenever the 
2fci?-CDW instability competes with the Coulomb driven three-fold charge order transition. The 
relevance of our results to the family of quarter-filled materials: S-(BEDT-TTF)2X is discussed. 



PACS numbers: 71.27.+a,74.40.Kb,71.30.+h,71.45.Lr 

I. INTRODUCTION 

Strongly correlated electron materials are often char- 
acterized by complex phase diagrams, reflecting an in- 
tricate interplay between magnetic, orbital, lattice and 
charge degrees of freedom. As these excitations cou- 
ple to the conduction electrons, the metallic state ex- 
pected in the absence of interactions has to compete with 
several ordered phases. Examples of these materials in- 
clude cuprate superconductors, nickelates, heavy fermion 
compounds, transition metal dichalchogenides, organic 
charge transfer salts and the iron based pnicitide super- 
conductors, all presenting various forms of magnetic, or- 
bital and charge order. Even when a metallic phase is sta- 
bilized, these systems are generally found to exhibit large 
effective mass enhancements and electrical resistivities 
violating the loffe-Regel-Mott (IRM) condition^ Sur- 
prisingly enough, such "bad" metallic behavior does not 
impede the emergence of superconductivity, but rather 
appears to be a prerequisite for the achievement of high 
critical temperature^. 

Charge ordered (CO) phases are commonly observed in 
the class of two-dimensional organic compounds 0-ET2X 
(ET=BEDT-TTF, bisethylenedithio-tetrafulvalene^pEl, 
and ascribed to the prominent role of electron-electron 
interactionJSI. At the non-interacting level, these com- 
pounds are predicted to be metals with 3/4- filled elec- 
tronic bands. The observation of electronic ordering im- 
plies that the magnitude of electron-electron interactions 
is comparable with the widths of the relevant electronic 



bands constructed from the tt molecular orbitals. In turn, 
the presence of such strong interactions raises questions 
about the nature of the metallic phase in these mate- 
rials, that should exhibit distinctive features of "corre- 
lated electron systems" in the Mott sense. The proximity 
to charge ordering instabilities, with the possible emer- 
gence of quantum critical points as the transition tem- 
perature is made to vanish, is also expected to strongly 
alter the physical properties of the metal. All these in- 
gredients should lead to me asura ble deviations from the 
usual Fermi liquid behavioi^SEIl^ in close analogy with 
heavy fermion system^ ^^-H^SJ. 

The minimal theoretical description of the electronic 
properties of 0-ET2X organic conductors is based on 
the two-dimensional extended Hubbard model (EHM) 
on a triangular lattice. Several theoretical studies have 
aimed at reproducing the different CO patterns real- 
ized in this class of materials, either within the frame- 
work of the EHM itself or its generalizations, including 
longer ranged electronic inte ractio ns and various types 
of electron-lattice interaction ^ I In the present work 
we focus on the following open issues: (i) how does the 
strength of the local Coulomb correlations modify the 
nature of the metallic phase as well as its CO instabil- 
ities, (ii) what are the effects of geometrical frustration 
in the electron motion arising from the triangular molec- 
ular arrangement, and (iii) how does the proximity to a 
given CO phase extend its infiuence onto the properties 
of the correlated metal, possibly leading to non-Fermi 
liquid behavior? 

In Sec. In] we set the minimal electronic model needed 
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FIG. 1: (Color online) (a) Arrangement of BEDT-TTF 
molecules in the conducting layers of 9-type ET crystals, with 
the corresponding transfer integrals and nearest-neighbor 
Coulomb interactions, (b) Threefold charge ordered phase 
in the triangular lattice considered in this work. The triangu- 
lar lattice is defined by the real space unit vectors ai = (1,0) 
defining the c direction and sl2 = (1/2, V3/2). 

for the study of electronic properties of 6'-(ET)2X com- 
pounds and provide a brief review of established theo- 
retical results. In Sec. |III| the model is solved by Lanc- 
zos diagonalization and the zero-temperature phase di- 
agram is obtained for different degrees of geometrical 
frustration. In Sec. |IV| the resulting metallic phases 
are explored through a finite-T Lanczos diagonalization 
calculatior|22ii24|_ theoretically explore the conse- 

quences of a quantum critical point (QCP) at a charge 
ordering transition driven by the quantum fluctuations 
associated with strong inter-site Coulomb repulsion. Our 
results are compared with a spinless calculation in order 
to assess the importance of the magnetic degrees of free- 
dom in the observed quantum criticality and make con- 
tact with the existing literature. The relevance of the 
present results to the physics of 0-type ET compounds is 
discussed in the conclusive Sec. |Vl 



II. MODEL AND METHOD 

1. Extended Hubbard model on the triangular lattice 

Quasi two-dimensional organic conductors (general 
formula A2B) are charge transfer compounds composed 
of alternating layers of conducting (donor) molecules A 
and insulating (acceptor) units B. They exhibit a large 
variety of molecular arrangements corresponding to dif- 
ferent polytypes classified by greek characters''^ The 
materials of the 6'-ET2X class have a triangular lattice 
structure, shown in Fig. [l] with an average of n = 3/2 
electrons per molecule, fixed by complete charge transfer 
between A and B units. Since this corresponds to a three- 
quarter filled electronic band, these materials should be 
normal metals in the absence of interactions. 

The electronic properties of 6'-ET2X materials are 



commonly described via the extended Hubbard Model 
(EHM), 

(ij)p(7 {ij)ciy 

+U'^ni^ni^ + Vp'^ niUj + V"c ^ n.,nj. (1) 

This model includes transfer integrals between the tt- 
orbitals of nearest-neighboring molecules in the conduct- 
ing plane, labeled by tp and according to the bond 
directions, a local (on-site) Coulomb repulsion energy U 
as well as non-local (nearest-neighbor) repulsion terms 
Vc,Vp (see Fig.ll|'^i. 

Early mean-field calculation^2llMl indicated that three 
types of striped patterns (vertical, diagonal and horizon- 
tal) are realized depending on the relative magnitude of 
the non-local Coulomb interaction parameters. These 
results were later confirmed by more advanced numeri- 
cal techniques that can properly account for electronic 
correlations, such as exact diagonalization (ED)^^'^ and 
density matrix renormalization group (DMRG)3il. 

More interesting from our persp ective is th e isotropi- 
cally interacting case, Vp ^ Vc = iMMimmm, There, 
because of the frustration of inter-molecular interactions 
induced by the triangular lattice geometry, an alterna- 
tive charge ordering pattern with three-fold periodicity is 
favored with respect to the (degenerate) striped arrange- 
ments, illustrated in Fig. [ijs. A more exotic situation is 
found in the limit of strong local Coulomb interactions 
(or, similarly, in a fully spin polarized electron system, 
i.e. for spinless electrons), where the constraint of no 
double occupancy on molecular sites converts this three- 
fold order i nto a parti ally ordered phase termed "pinball 
liquid" (PL )'^ ' ^' '33135i. ^^lis state shows a three-sublattice 
structure with the same symmetry as the threefold phase, 
in which the carriers of one sublattice are essentially lo- 
calized as a Wigner-crystal (pin), with the remaining 
charges (balls) forming an itinerant liquid on the intersti- 
tials. It is not clear at present how the transition between 
these two qualitatively different forms of threefold order 
takes place as a function of the local Coulomb repulsion 
U. This issue will therefore be thoroughly discussed here. 

In addition to the effects of the local electronic cor- 
relations, we are interested in the effects of geometrical 
frustration in the electronic motion, that are triggered 
by the strongly directional 7r-overlaps between neigh- 
boring moleculeJ^. This issue is of particular impor- 
tance to actual materials, as the relative values of the 
transfer integrals tc and tp can be tuned experimen- 
tally by applying pressure or by chemical substitution, 
which m odifie s the relative angles between neighboring 
molecule^SlllZl^ general observation, negative val- 

ues of the ra tio te/ ^D produce the highest charge ordering 
temperature while vanishing or positive values lead 
to glassy (X=CsCo(SCN)4, X^CsZn( SCN)4l^ or even 
superconducting (X=l3) ground state^^2ll6l39]_ pj-om a 
more theoretical point of view, how the system evolves 
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from a perfectly isotropic triangular lattice at tc = ±tp 
to a square lattice at = remains an open issue. 

Transfer integrals obtained through the Hiickel approx- 
imation in quarter-filled 0-ET crystals are in the range: 
-0.5 < tc/tp < 1.5, with tp « -0.05 - 0.1 e\^. These 
values generally differ from the ones extracted from op- 
tical reflectivity and de Haas Van Alphen experiments^Sl 
for each specific crystal. On the other hand, Coulomb 
repulsion energies in organic molecular crystal^^ have 
been estimated by calculating the screening corrections 
to the bare repulsion energies of the isolated molecules, 
Uq and Vo, obtained from ab initio calculations^^. These 
calculations lead to Hubbard parameters of the model 



and Vp^Vc^ U/2, 



Eq. ([T]): [/ ~ Uo/2 ~ (15 - 20)\t 
with a bandwidth W ~ (8 — 9)|tp|. These Coul omb ener- 
gies are larger than assumed in previous work^^SllIIlIIlIl 
[/ ~ (8 — 10)|tp|, and V ^ {1 — 3) |tp | , as extracted from 
optical reflectivity measurementa^. The degree of uncer- 
tainity in the microscopic parameters implies that a gen- 
eral understanding of the model and its phase diagram 
in the full parameter space U, V and tc/tp is essential. 
This is the main focus of the present work. 



2. Fmite-T Lanczos approach 

We perform ED calculations through a flnite-T Lanc- 
zos algorithm with periodic boundary condition^22i23j^ 
The large number of excited states inherent to the many- 
body problem which are needed to evaluate statistical 
sums is cutoff by keeping only a small number of low lying 
states at each temperature. This is performed through an 
Arnoldi algorithnpS which reduces the size of the Hilbert 
space enormously. The accuracy of the method is re- 
stricted to temperatures which are not too low, i.e. not 
lower than the energy of the lowest excitation of the quan- 
tum many-body system. For the method to be practical 
T should not be too large that one needs to keep too 
many states in the statistical sums. Finite-size effects are 
somewhat reduced by the effect of temperature and the 
method is quite reliable for extracting integrated proper- 
ties. Instead, spectral properties such as optical conduc- 
tivity and photoemission spectra are prone to large finite 
size effects, and will not be analyzed here. 

Due to the high computational demand of the finite-T 
algorithm, the calculations are performed on an Ng = 
12 site cluster. In principle, a larger Ng = 18 cluster 
whose geometry is also suitable for reproducing the three- 
fold CO pattern could be used at T = 0. However, we 
shall not consider such case because a 3/4-filling implies 
a different number of spin up and spin down electrons, 
while the ground state is expected to be in a S* = 5*2 = 
state. 

We characterize the physical properties of the different 
phases based on the following quantities, that are acces- 
sible through finite-T Lanczos calculations. 

(i) Charge correlation function: The charge structure 
factor signalling the possible occurrence of a charge or- 



dered state in the system is evaluated at finite-T through 

1^ 



^ninj\m). (2) 



Here Z — g-0Em jg partition function of the 
system and /3 = l/fc^T. A charge ordered state with 
modulation Q is signalled if C(Q) is finite in the ther- 
modynamic limit. The three- fold ordering corresponds 
to a charge density modulation with wavevector Q — 
(27r/3, 27r/\/3), which lies at the corner of the hexago- 
nal Brillouin zone, see Fig. 11 (all corners are equiv- 
alent, being connected by reciprocal lattice vectors or 
time-reversal symmetry). An accurate numerical deter- 
mination of the phase boundaries should rely on a proper 
finite-size scaling of the results. While this is prohibitive 
for the fermionic system under study due to the rapidly 
increasing size of the Hilbert space, the ordering transi- 
tions can still be identified as the locus of steepest varia- 
tion of charge correlations upon varying the microscopic 
parameters of the model. 

(ii) Kinetic energy: This quantity provides direct in- 
formation on how the motion of the charge carriers is 
slowed down by interactions. It can be evaluated with 
high accuracy from the flnite-T Lanczos diagonalization, 
because it results from a quantum mechanical and ther- 
mal average over a huge number of states. By normal- 
izing it to a reference non-interacting value Kq^ it gives 
valuable information on the degr ee of electronic corre- 
lations in the many-body systerrp2144|_ Under suitable 
assumptions, this quantity can be compared with optical 
absorption experiments in actual materials via the f-sum 
rule^^ 

The kinetic energy is evaluated from the following ther- 
mal average. 
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where \m) is the total set of eigenstates of the system 
with energies E„-^. 

(iii) Double occupancy: It is useful to analyze the num- 
ber of double occupancy per site in the lattice which 
reads: 



(4) 



and is different for the different phases analyzed. For 
example, it is a key quantity in the analysis of the Mott 
transition in the half-filled Hubbard model since d is sup- 
pressed in the Mott insulator, which allows to determine 
the critical Coulomb coupling. In model Eq. ([l]), c? is 
helpful for characterizing the different possible CO states 
for different U and V . 

(iv) Specific heat: From the total energy of the system, 
E = (if), we can obtain the specific heat by taking the 
derivative with respect to the temperature, T: 
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FIG. 2: (Color online) Phase diagram obtained at T = 
from numerical diagonalization oi a. Ns = 12 cluster. The dif- 
ferent shaded areas correspond to the pinball (PL) phase for 
tc/tp — 1, 0.5, 0, —0.5 (squares, circles, triangles, diamonds re- 
spectively, from light to dark blue). The case tc/tp = — 1 has 
no pinball phase in the explored U range. The red line cor- 
responds to a direct transition from the homogeneous metal 
(HM) or pinball liquid (PL) to the threefold charge ordered 
(3C0) state. 

Unless otherwise specified, we use units such that 
kg = h = 1. The finite temperature method recovers 
the ground state properties by taking the hmit: /3 oo. 
In practice this is achieved for /3 = 50 — 100 for the var- 
ious U and V explored across the whole phase diagram. 
Typically about 30 to 50 terms are kept in the evaluation 
of the statistical sums over the excited states |m) with 
corresponding energies Em- 



III. PHASE DIAGRAM AT T=0 

The zero temperature phase diagram of the model Eq. 
([T]) in the {U, V) plane is shown in Fig. [5] The phase 
transition lines are determined using three alternative 
methods, that all give coincident results: (i) via the evo- 
lution of the charge correlation function C(Q) calculated 
at the three-fold wavevector. Fig. |3] (ii) by tracking di- 
rectly the charge ordering patterns that develop in real 
space, Fig. [4j and (iii) by analyzing the fidelity between 
groundstates at different values of the microscopic pa- 
rameters. Fig. [5] as introduced below. 

We are interested here in the charge ordering instabil- 
ities driven by the inter-site repulsion V. Our numerical 
results confirm the existence of three distinct phases: a 
homogeneous metal (HM) at low V, a three-fold charge 
ordered phase (SCO) at large V, and an intermediate 
"pinball liquid" (PL) phase emerging at large values of 
U. The most striking effect in Fig. [2] is that the region 
of the homogeneous metallic phase is strongly reduced 
upon increasing the tc/tp ratio, evidently due to a cor- 
responding stabilization of the competing pinball liquid 



phase. We note that the homogeneous metal is always 
the ground state at = independently of the strength 
of the local repulsion U. This can be rationalized by the 
fact that in the absence of nearest-neighbor interactions, 
at = 3/2 the holes {uh = 2 — n=l/2) can effectively 
avoid each other when moving along the lattice. 

A. Characterization of the different phases 

Fig. [3] reports the evolution of the charge correlation 
function, the kinetic energy and double occupancy as a 
function of V, along vertical cuts in the phase diagram 
corresponding to U/tp — 5 and U/tp — 40. Different 
curves correspond to different values of tc/tp (upper and 
lower panels, respectively). 

1. Small U: HM to SCO transition 

At low J7, the instability towards the threefold charge 
ordered phase is signaled by a sharp jump in the corre- 
lation function (Fig. [3^), starting from a small constant 
value in the homogeneous metal. The locus of the 3C0 
transition shows an appreciable dependence on geometri- 
cal frustration: the homogeneous metal is rapidly desta- 
bilized for positive values of tc/tp^ in marked contrast 
with the weaker (and opposite) variations expected from 
an RPA analysis valid in the weakly correlated limif^ 
(see Appendix [A|. From the phase diagram of Fig. [2] it 
is quite clear that this trend is governed by a mechanism 
that extends from the strongly correlated limit U ^ tp 
down to the lowest values of U. The emergence of an 
intermediate plateau in the charge correlation function, 
clearly visible in the data at tc = tp in Fig. [3^, is also 
reminiscent of the situation encountered at U/tp = 40 
(see below). These observations suggest that the pres- 
ence of geometrical frustration, tc/tp > 0, strongly en- 
hances the role of electronic correlations. The pinball 
phase characteristic of strong U is stabilized at tc = tp 
despite a relatively low nominal value U/tp — 5. 

We note that the nature of the ordering transition 
changes in the opposite limiting case tc — —tp, where 
the sharp jump in the correlation function is replaced by 
a smoother evolution, possibly due to the competition 
with an incipient nesting instability (Appendix [A| . 

The behavior of the charge correlation function is di- 
rectly mirrored in the other physical quantities shown 
in Figs. [3j The kinetic energy (Fig. [3]:) jumps at the 
phase transition from an essentially free-electron value, 
K/Kq > 0.9, to a value that is reduced by the open- 
ing of the charge ordering gap. At the same time, the 
double occupancy (Fig. [sjs) undergoes a marked increase 
towards the value d — 0.66 of the fully formed SCO: the 
charge is ordered into three sublattices with average oc- 
cupations UA = riB = 2 and nc = 1/2 (see Figjij large 
V region), so that each of the two charge rich sublattices 
contributes dA — ds = 1/3 to the average double occu- 
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FIG. 3: (Color online) (a)-(b) Charge correlation function for the threefold wavevector Q = (27r/3, 27r/\/3) as a function of 
V/tp (a) in the weakly correlated U — 5tp and (b) the strongly correlated limit U — 40tp. (c)-(d) Average normalized kinetic 
energy K/Kq as a function of V/tp. (e)-(f) Double occupancy probability. The plots in (b) have a vertical offset of ±0.02 for 
tc — ±tp and of ±0.01 for tc ± 0.5tp for clarity. 



pancy. The fact that the double occupancy in the ho- 
mogeneous metal is suppressed from the non-interacting 
value d = (n/2)2 = (3/4)2 ^ 0.5625 indicates the pres- 
ence of moderate electronic correlations. 



2. Large U : HM to PL transition 

A richer situation is found in the large U regime. First 
of all, the homogeneous liquid is characterized by a total 
suppression of double occupancy: introducing the double 
occupancy of holes dh — ((1 — ni^){l — ni\^)) = 1 — n ± c? 
and setting d/j = we obtain d — n — 1 = 1/2, which is 
actually observed in Fig. |3f. Furthermore, the presence 
of both local and non-local Coulomb interactions hinders 
the particle motion, resulting in a marked reduction of 
the kinetic energy upon increasing V . An approximate 
expression for its ^-dependence is: 

k'^"'') ^{I^AV^)k\j"^"\ (6) 

where K^^^'' is the value atV — 0. This dependence 
is consistent with a previous slave-bosorPS! calculation of 
the metallic phase formed by spinless particles on a d- 
dimensional hypercubic lattice, which is compatible with 
the data of Fig. |3jl. Because the ordering instability is 
pushed to large values of V for negative values of tc/tp (cf. 



Fig. [2]), the homogeneous metal that is so revealed can 
acquire quite a strongly correlated character, as testified 
by a kinetic energy ratio that decreases down to K/Kq — 
0.5 before the onset of charge order. 

The numerical data of Figs. [SJd and d show quite 
clearly that an intermediate phase emerges between the 
homogeneous metal and the SCO, that we associate 
with the pinball liquid phase introduced by Hotta and 
coworkera22E3EZl xhe PL is a partially ordered phase 
with a three-sublattice structure, in which the carri- 
ers of one sublattice (pins) are localized as a Wigner- 
crystal and the remainder (balls) form an essentially non- 
interacting liquid within the resulting hexagonal lattice 
(Fig.lb). 

To get a further insight about the different broken sym- 
metry phases, we have calculated the static density pro- 
file in the presence of a local perturbation breaking all 
the translations of the lattice, but respecting the 7r/3 
rotation. In this way, instead of obtaining a uniform lin- 
ear combination of all symmetry related crystal states, 
the system selects one crystal state favored by the per- 
turbation, giving access to a real-space snapshot of the 
broken-symmetry ground state. The basis of the method 
employed for getting the real space snapshot of the lo- 
cal densities is the following: we add a local potential on 
certain sites related by rotations in such a way that only 
translations are broken. The value of the local defects is 
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FIG. 4: (Color online) Average electron densities in the elec- 
tron rich (hole poor, ua = ub > n, upper curves) and electron 
poor (hole rich, nc < n, lower curves) sublattices, as obtained 
from Lanczos diagonalization at U/tp =40 in the presence of 
a weak translational symmetry breaking potential. We take 
5 — +0.05tp on the sites of the hole rich sublattice in or- 
der to select one realization of the threefold symmetry. Color 
codes are the same as in Fig. [2] (from light to dark blue: 
tc/tp = 1,0.5,0,-0.5, red: t^/tp = -1). 



few percents of the hopping term. In the Hamiltonian, 
we simply add the perturbation term: SJ^i ("-it + "-4)- 
This method has previously been used for distinguishing 
the exact nature of two phases breaking the translational 
and the rotational symmetries differently *^. Since the ad- 
ditional term is kept very small, it only corresponds to 
a perturbation and do not change the main property of 
the ground state. In this paper, we do want to keep the 
rotational symmetry safe, we hence chose 4 sites, in the 
12-site cluster, in such a way that the system is still ro- 
tationally invariant. The perturbation has to be lower 
than the typical energy scale and we put S — +0.05tp 
on the sites of the hole rich sublattice. Results are de- 
picted in Fig. |4j Expressing the charge densities in terms 
of holes, Uh = 2 — n = 1/2 and starting from the SCO 
phase {fih^A = nh^B = and nh,c = 3/2, Fig. [I]), by pro- 
gressively increasing the local repulsion ?7, it becomes 
energetically unfavorable to acommodate more than one 
hole per molecule. Part of the hole density will then 
tend to spill out of the hole-rich sites in order to prevent 
double occupancy, resulting in rih^c = 1- The average 
charge density in the three sublattices, ua — tlb — 7/4, 
nc = 1, illustrated in Fig. |4] indeed corresponds to a sce- 
nario where one sublattice is occupied by localized holes 
{nh,c = 1 leading to nc = 2 — n^.c = 1) with the re- 
maining particles equally spread in the interstitial sites 
{nh,A = nh.B = iT-baii = 1/4 hole per site, leading to 
nA — nB = 2 — nbaii = 7/4). The critical value for the 
transition from SCO to PL can be readily estimated to 
be J7 ~ SV^ from electrostatic considerations alone (see 
Appendix [B]) , and is therefore independent of tc/tp. 

Notably, in the PL phase the physical quantities whose 
V dependence is depicted in Fig. [3] form well defined 
plateaus, suggesting that they are locked as a conse- 



quence of the spontaneous separation between localized 
and itinerant charges. For example, the charge correla- 
tion function tends to the value C(Q) = n^/3 ~ 0.08 
instead of the full C(Q) = n^ obtained at complete or- 
dering in the large V limit, corresponding to the fact 
that only 1/3 of the particles participate to the ordering 
phenomenon. 

From our numerical results in Fig. [3|1, the kinetic en- 
ergy of both the homogeneous metal a,t V < VpL and the 
SCO phase aXV > V^co is found to depend only weakly 
on the degree of geometrical frustration tc/tp. On the 
other hand, the kinetic energy forms a plateau within 
the PL phase, at a value which is strongly dependent 
on tc/tp: the absolute value of the (negative) kinetic en- 
ergy in the PL phase is lowest at — —O.htp and its 
magnitude steadily increases as the frustration ratio is 
increased to tc = tp. Our data therefore suggest that it 
is this gain in kinetic energy for positive values of the 
geometrical frustration ratio that is responsible for the 
stabilization of the pinball liquid phase against the ho- 
mogeneous metal observed in Fig. [2p2I We note that a 
kinetic energy driven mechanism for the PL transition 
was also pointed out in Ref. 34, Since the effective fill- 
ing associated with the itinerant balls on the hexagonal 
lattice is only njyaii — 1/8, Coulomb interaction effects 
are small. Indeed, we have checked that the kinetic en- 
ergy of the balls coincides with the kinetic energy of the 
corresponding non-interacting tight-binding model on a 
hexagonal lattice, at filling n = ntaih as a function of the 
ratio tc/tp. 

B. Fidelity analysis 

In order to characterize the quantum phase transi- 
tions of our system, the simple yet powerful concept of 
fidelity^ is considered here (Fig. [sj . Initially introduced 
in quantum information, the fidelity has revealed suc- 
cessful in determining superfluid-insulator transitions of 
the Hubbard modeP^JM] -p-j^^ idesu behind the fidelity is 
very simple; it consists of computing overlaps of ground 
states (GS) |-0o) at different values of the microscopic 
parameters. At a QCP, even the smallest change of the 
parameters can have dramatic effects in some of the ob- 
servables. This is encoded in the GS properties, hence, 
the overlap is expected to strongly react and indicate the 
locations of the phase transitions. We define the fidelity 
F as: 

Fvy,^\{MV')\MV))\. (7) 

Obviously, for V = V , the fidelity should be Fyy = 1- 
Typical results are depicted in Fig. [5] at [/ = 40ip, and 
for five values of t^tp = ±1.0, ±0.5 and 0.0. 

It is surprising to see how the fidelity is indeed able 
to pinpoint the phase transitions. For each of the values 
of tc/tp, we exactly recover the transitions obtained by 
more standard methods in the preceding Sections. More 
information is available, however. First, there is always 
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FIG. 5: (Color online). Fidelity Fv,v' = \{Tpo{V')\ipo{V))\ for various values of tc/tp, from left to right 1.0,0.5,0.0,-0.5 
and —1.0 and U = 40|tp|. There is a perfect agreement with the phase transitions obtained by standard quantities (order 
paramaters, kinetic energies, suceptibilities). The case of zero overlap is represented in black, while the blue- violet (diagonal 
squares) represents overlap 1 {V^V'). The non-uniform colored zones stand for fluctuating ground states possibly related to 
a strong degeneracy. 



some region where F is strongly fluctuating, even though 
weU dehmited in the (V, V) plane (non uniform colored 
zones). These fluctuations can be due to a large degen- 
eracy of the GS hence corresponding to a same order but 
with destructive interferences. The second important in- 
formation is visible in the case tc/tp — 0.5, where the 
fidelity indicates the existence of two distinct phases in 
the pinbaU region [namely 5.20(2) < V/tp < 11.00(2) 
and 11.00(2) < V/tp < 13.00(2)], suggesting a possible 
ordering of the mobile charges. In fact, a slight change 
in the kinetic energy (Fig. [sji) appears at this transi- 
tion, but the other quantities seem to be insensitive to 
it. The fact that the transition within the pinball phase 
is not detected by the charge correlation function (Fig. 
[Sja) nor by the average densities (Fig. |4]) tends to show 
that their properties remain extremely close. Neverthe- 
less, the fidelity F allows a precise determination of the 
sub-phases. 

C. Spinless model 

We now consider the spinless version of the model 
Eq. ([T]), which has been discussed extensively in the 
literatur^^mmZIlll, Since the spinless model only con- 
tains charge degrees of freedom, by comparing it with 
the spinful model we can obtain useful information on 
the relative role played by charge fiuctuations as com- 
pared to the spin fluctuations. 

The spinless model reads: 

H = tpY^ {h\hj + iJ.c.) + 51 + ^-c-) 

+ Vp'^ fiiTij + Vc ^ niTij, (8) 

{ij)p (ij>c 

where we have changed the sign of the hopping integrals 
to deal explicitly with holes. As in the preceding Sections 
we consider the case Vc — Vp = V for different values 
of the tc/tp ratio. Importantly, for spinless particles the 
physical situation of one hole per two sites implies a half- 
filled band, which gives rise to a spurious particle-hole 




V/t 

p 



FIG. 6: (Color online) Dependence of the kinetic energy 
on V for the spinless and spinful extended Hubbard models. 
Different tc/tp ratios are compared in the two situations. As 
expected, the properties of the spinless model are independent 
of the sign of tc/tp, in contrast to the spinful model. 



invariance that is absent in the spinful case at 3/4-filling. 
Therefore, the thermodynamic properties as well as the 
phase transition lines become invariant under a change 
of sign of tc/tp. 

In Fig. [6] we show the ^-dependence of the kinetic 
energy for different tc/tp ratios compared with the spinful 
model. As expected the figure shows how the kinetic 
energy in the spinless model does not depend on the sign 
of tc/tp (black curves). The spinless model is able to 
recover the qualitative trends of the full model for tc/tp > 
0: the critical values at which the transition from the HM 
to the PL occurs are respectively Vp^ ~ 4.5ip and 8.5tp 
for tc/tp ~ 0.5 and in the spinless case, to be compared 
with Yph — 5.5tp and 7.5tp in the spinful case at U — 
AOtp (Fig. |6]). However, because of the artificial particle- 
hole symmetry, for tc/tp < the locus of the transition 
is completely inconsistent with the spinful case, and the 
stability of the PL is widely overestimated. 

We see from Fig. [6] that the effect of non-local inter- 



8 



actions V on the renormalization of the kinetic energy 
in the homogeneous metal is very similar in the spin- 
ful model at large U and in the spinless model, both 
being compatible with the quadratic F-dependence dis- 
cussed above. We therefore conclude that the different 
behaviors observed in the two models at tc/tp > are 
a direct consequence of the different kinetic energies at 
the non-interacting level, that results from the spurious 
particle-hole symmetry acquired by the spinless version. 

The results presented here indicate that the charge 
rather than spin correlations dominate the renormaliza- 
tion effects on metallic properties approaching the charge 
order transition. However the spin multiplicity enters 
(indirectly) via the geometrical frustration, that is not 
treated correctly in the spinless model. A realistic spin- 
ful calculation therefore appears to be necessary to prop- 
erly address the physics of 6'-(ET)2X salts, where tc/tp is 
a key parameter in determining the experimental phase 
diagranPl2i. 

IV. CORRELATED METAL AT FINITE 
TEMPERATURES 

Here we analyze the properties of the homogeneous 
metallic phase at finite temperatures close to the QCP. 
A temperature scale emerges, that we denote T*, above 
which the kinetic energy departs from Fermi liquid be- 
havior and the specific heat coefficient goes through a 
maximum. We interpret T* as a renormalized Fermi tem- 
perature, that generally drops to zero at the approach of 
the QCP. Such behavior is typically found for geomet- 
rical frustration tc ^ ±ip. On the contrary, in cases 
in which there are competing Fermi surface instabilities, 
especially in the perfectly nested case tc = ~tp, the T* 
phenomenon is much weaker, and hardly affects the prop- 
erties of the electron liquid. The data in that case are 
compatible with a Fermi temperature that remains finite 
right close to the QCP, possibly indicating that a first 
order transition may be occurring. 

In the following paragraphs we focus specifically on 
the realistic value U — 15ip, but the qualitative features 
presented here are unchanged for different large values of 
U. 



A. Non-Fermi liquid behavior close to CO 

1. Kinetic energy 

The temperature dependence of the average kinetic en- 
ergy K normalized to the non-interacting value Kq [U = 
y = 0) is shown in Fig. [t] for tc/tp = 1,0.5,0,-0.5,-1. 
The different curves in each panel correspond to differ- 
ent values of the intersite Coulomb repulsion, V , across 
the charge ordering transitions. The ratios tc/tp — 
0.5, 0, —0.5 (Fig. [7jD, c and d) all show marked depar- 
tures from the quadratic temperature dependence K — 
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FIG. 7: Temperature dependence of the total kinetic energy 
at [/ = 15tp normalized to the non- interacting value Kq, for 
tc/tp = -1,-0.5,0,0.5 and 1. Arrows correspond to inflec- 
tion points of the curves, indicating the temperature scale T* 
corresponding to the breakdown of the quasiparticles. Open 
circles correspond to the HM phase, diamonds to the PL (full 
lines) and squares (dashed lines) to the SCO phase. 

Kt=o — BT^ characteristic of conventional metals, oc- 
curring in the HM phase above a certain temperature 
(the approximate locus of the inflection points is indi- 
cated by arrows). We denote it as T* and take it as an 
estimate of the renormalized Fermi temperature, govern- 
ing a crossover to non-Fermi liquid behavior. Clearly, 
T* is progressively reduced upon approaching the charge 
ordering transition and vanishes at the critical point. 
In this respect, our data in the experimentally relevant 
case U = 15tp do not show qualitative differences be- 
tween the transition to the threefold charge order ob- 
tained for tc/tp — — 0.5;0 and that to the pinball phase 
for tc/tp = 0.5 (cf. Fig. [2]): in both cases the temper- 
ature scale T* appears tooe entirely controlled by the 
approach to the zero-temperature ordering transition, in- 
dicating the possibility of quantum critical behavior at 
finite temperatures around the zero-temperature phase 
transition. 

In the cases in which tc = ±tp, the emergence of a 
temperature scale T* is much less clear [Fig. [t] (a) and 
(e)] as the temperature dependence of the kinetic energy 
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is smooth within the whole homogeneous metalhc phase 
except very close to the transition. Our RPA analysis 
presented in Appendix |X] shows that in these cases there 
is a competing CDW instability, that could indeed be 
masking the quantum critical behavior associated with 



the SCO transition. For tc/tp 



1 the CDW is driven 



by the perfect nesting of the Fermi surface whereas for 
tc/tp = 1 a mixed CDW/CO phase induced by both nest- 
ing tendencies and strong Coulomb repulsion coexist. 

The perfect nesting of the Fermi surface for tc/tp = 
— 1 occurs at the wavevectors Qi? — (±7r, ±7r/\/3) and 
(0, ±27r/-\/3). Such nesting instability is dominant at 
weak [/, where it results in a striped charge modulation 
that dominates over the threefold charge order discussed 
above (this striped order is analogous to the checkerboard 
pattern obtained in the square lattice). Our data sug- 
gests that even in the presence of a sizable local Coulomb 
repulsion, U — 15tp, that prevents the stabilization of 
such stripe order, an incipient nesting instability is strong 
enough to destroy the quantum criticality around the 
QCP. This conclusion is based on the observation that 
there is no clear signature of the vanishing low tempera- 
ture scale, r*, at the QCP and there is no clear evidence 
of the 'bad' metallic behavior found for other tc/tp ratios. 
However, there is a weak T-dependence of the kinetic en- 
ergy that vanishes at a critical value Vc — Qtp as can 
be observed from the data of Fig. [t] (e). Such critical 
value, Vc, is found to be consistent with the critical value 
obtained from the SCO charge correlations calculated be- 
low. 



2. Charge correlations 
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FIG. 8; Charge correlation function versus the temperature 
for different V/t in the CO phase a.t U = 15t for tc/tp = 
1,0.5,0,-0.5,-1 Arrows correspond to inflection points of 
the curves, indicating the ordering temperature Tco- 



The emerging QCP scenario can be further appreci- 
ated by studying the evolution of the charge ordering 
transition Tco vs. temperature. This can be obtained 
by tracking the steepest variation of the charge corre- 
lation function, C(Q) (arrows in Fig. [8^-d). In Fig. [o] 
we report Tco together with the T* extracted from fig- 
ures [T] showing a common behavior in proximity to the 
CO instability. As stated in the preceding paragraph, 
the case tc/tp = — 1 exhibits a different behavior, with 
no visible T* approaching the QCP. The charge correla- 
tions (Fig. [8^) also exhibit a qualitatively different be- 
havior in the perfectly nested case tc/tp — —1, with a 
mild non-monotonic temperature dependence showing a 
maximum at intermediate temperatures which suggests a 
'reentrant' ordering transition. For this case, there is no 
clear indication of 'bad' metallic behavior in the kinetic 
energy and there is clear evidence of 'reentrant' behav- 
ior in the SCO transition. This 'reentrant' behavior dis- 
appears at around the critical value Vc — Qtp as can be 
noted in Fig. |8|i. It can be noted that a slightly reentrant 
behavior can be extracted from the data at tc/tp = —0.5. 
The presence of a reentrant behavior at negative values 
of tc/tp is confirmed by our RPA analysis, and is strongly 
reminiscent of what is commonly observed in the EHM 



on the square lattice (cf. Fig. 1 in Ref.^''), which also in 
that case is ascribed to the competition of the CO phase 
with a Fermi surface nesting instability. We have added 
to the phase diagram for the case tc = —tp the 2kp- 
CDW instability and may also be present in other cases. 
However, the limited wavevector resolution of our small 
cluster calculation does not permit an accurate determi- 
nation of the stability of the 2/ci?-CDW phase. In fact, for 
tc — —O.btp, we may also expect that CDW instabilities 
occur in the proximity of the SCO instability. 



3. Specific tieat 

Further insight on the anomalous properties of the ho- 
mogeneous phase can be gained by exploring thermody- 
namic properties such as the temperature dependence of 
the specific heat coefficient Cy/T on approaching the 
QCP. In a Fermi liquid at low T, this quantity measures 
the effective mass enhancement of the quasiparticles. In 
we compare Cy/T for tc/tp = 0.5,-0.5,-1. In 
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Fig. _ 

this way we compare the behavior of the specific heat of 



a system across the SCO {tc 



-0.5t„) with a system 
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FIG. 9; (Color online) Color plots summarizing the QCP behavior a.t U = 15fp for the representative cases tc/tp = 0.5, —0.5 
and —1. The dashed curve is the renormalized Fermi temperature T* extracted from Fig. [t] The full line is the charge ordering 
temperature Tco extracted from Fig. |8] The color gradients (blue and red respectively) are derived from the kinetic energy 
and charge correlation data of Figs. [7] and [5] 
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FIG. 10: (Color online) Temperature dependence of effective mass enhancement for different tc/tp ratios for U = 15tp. A 
different behavior in the T-dependence is found for the case tc — —tp for which T* (indicated by the vertical arrows) which 
indicates a departure from the quantum critical behavior found in the other non-nested Fermi surface situations. 



across the PL transition {tc/tp — 0.5). For completeness 
we also analyze the t^/tp = — 1 case in which perfect 
nesting exists which can be compared to the other two 
cases. For both tc/tp = ±0.5 we find that a peak in 
Cv/T develops at T* with T* — >• on approaching the 
QCP asV ^Vc. This indicates that both the drop of T* 
and the effective mass enhancement occurring in proxim- 
ity to the QCP are consistent with the phase diagram of 
Fig. 9. Both effects are key signatures of the presence of 
a QCP together with the 'bad' metallic behavior arising 
around it. 

The case tc — —tp deserves special attention. For 
Coulomb repulsion energies up to about V = bAtp, there 
is a moderate increase of the effective mass enhancements 
and a moderate shift of the peak to lower temperatures 
in contrast to previous cases. However, increasing V fur- 
ther leads to different behavior with a rapid increase of 
the effective mass enhancement and a shift of the peak to 
zero which indicates the proximity to a SCO transition 
at Vc — Qtp. The shift of the peak in the specific heat 
with V is plotted in the phase diagram of Fig. |9] together 
with the transition line to the SCO. The two lines merge 



and a clear 'reentrant' behavior of the SCO transition is 
observed which is ascribed to the presence of the compet- 
ing Fermi surface nesting instability. The origin of the 
specific heat peak is unclear since it does not separate the 
HM from the 'bad' metal and from Fig. [9] it is clear that 
the peak position has a different qualitative behavior for 
tc/tp — —1 than for other tc/tp ratios. In fact, there is 
no 'bad' metallic behavior in the HM phase as discussed 
previously. However, the fact that the effective mass is 
enhanced may indicate a transition to the 2fc/-CDW. 



B. Connection with QCP physics 

It is worth analyzing our numerical results from 
the pe rspec tive of the standard Moriya-Hertz-Millis 
(MHMJ^^ theory of quantum critical points (QCP). In 
principle, this theory could be appropriate to the tran- 
sition from the HM to the SCO phases considered here 
since both are metallic. However, one should keep in 
mind the limitations of the MHM theory. First, it is 
based on a weak coupling perturbative expansion around 



11 



the QCP. Second, the MHM theory is not apphcable for 
ordering transitions which are driven by "2k/" Fermi 
surface instabihties. Hence, the MHM theory is ap- 
phcable only to systems in which there is no nesting 
at all and for which the ordering wavevector Q is not 
an extremal vector of the Fermi surface, i. e. loosely 
speaking Q ^ 2k f. Our model satisfies the latter con- 
dition whenever tc/tp ^ ±1 for the SCO wavevector: 
Q = (27r/3,27r/y3) ( see Fermi surfaces and discussion 
in Appendix|A]). In these cases, the two-dimensional ver- 
sion of MHM, d = 2, with a dynamical scaling dimension, 
z = 2 (corresponding to nearly antiferromagnetic metals 
with effective dimension at the QCP: D = d + z = 4) is 
relevant. This situation corresponds to a marginal case 
which contains dangerously irrelevant operators in the 
renormalization group sense which can destroy the hy- 
perscaling at the QCP. In this marginal case, d = z = 2, 
the Fermi liquid phase is bounded in the phase diagram 
by the condition: T < r, where r — \V — Vc\ quan- 
tifies the proximity to the QCP from the metallic side 
of the transition and Vc is the critical value at which 
CO occurs. The quantum critical region is bounded by 
the condition T > r obtained from the scaling behavior 
of the renormalization group (RG) equations. A simi- 
lar linear dependence with r is found in the boundary of 
the ordered phase which displays a critical non-Gaussian 
behavior around it due to the Coulomb interaction. As 
summarized in Fig. [9j for tc/tp ^ ±1 the suppression 
of T* at the QCP and the 'bad' metallic behavior ob- 
tained at finite-T from our numerical calculations are 
qualitatively consistent with the MHM predictions for 
the marginal d = z = 2 case. 

On the other hand, the extremal values: tc/tp — ±1 
deserve special consideration. In the particular case: 
tc = —tp, the system has perfect nesting at the ordering 
wavevector: Qp = (tt, 7r/-\/3), (see Fig. 11 in Appendix 
[A| ) which describes diagonal stripe order in real space. 
A RPA analysis on the model shows that a CDW insta- 
bility at Qp exists at small but finite U and V, which 
competes with the SCO with Q = {2n / 3 , 2n / ^/3) . Such 
coexistence/competition between CDW and CO instabil- 
ities has also been found in RPA studies of the extended 
Hubbard model on the square lattice-Sl. The situation in 
which two instabilities coexist — a nesting driven CDW 
with ordering vector Qp and a Coulomb driven instabil- 
ity with ordering vector Q — has not been addressed in 
general at the level of the MHM approach. Our anal- 
ysis shows that at moderate values of C/ > 7tp nesting 
instabilities are washed away and the charge ordering 
transition is Coulomb driven. However, the competition 
washes out the T* phenomenon and we find no clear ev- 
idence of 'bad' metallic behavior around the QCP in the 
case tc = —tp. This coincides with the breakdown of the 
MHM approach when nesting is present in the lattice. 

We may speculate, based on our numerical analysis, 
that in the perfectly nested situations quantum critical- 
ity is destroyed and a first order transition occurs. In- 
deed, a somewhat related renormalization group (RG) 



approach!^ to 2fci?-density wave quantum phase transi- 
tions in which curved Fermi surfaces with parallel tan- 
gents at two points of the Fermi surface connected by 
2kp are considered has found that critical fluctuations 
strongly influence the fermions on the Fermi surface and 
that the feedback effect of these fluctuations can destroy 
the second order quantum critical point turning it into a 
first order transition. Only in the special case in which 
Q = G/2, with G being a reciprocal lattice vector, a 
second order quantum phase transition is recovered. 

In actual quarter-filled organic materials, 0-ET2X, for 
which hopping ratios tc/tp ^ — 1 the MHM theory may 
be relevant. Many of the predictions for thermodynamic 
and transport properties in the quantum critical regime 
above the zero temperature QCP could be then experi- 
mentally checked. One important prediction of the MHM 
theory for d = 2: = 2 is the anomalous temperature de- 
pendence of the specific heat : 



C 



Tlnl/r, T<r (9) 
Tlnl/T, T>r, (10) 



where r describes the proximity to the QCP. On the other 
hand, in clean nearly charge ordered two-dimensional 
metals the resistivi ty ar ound the 'hot' spots in the quan- 
tum critical regimdi^lSZI reads: 



p (X T\ T^r 
T, T>r. 



(11) 
(12) 



However, the resistivity is shortcircuitcd by the contri- 
bution of electrons at the 'cold' spots since the scat- 
tering is small around these partJ^. Hence, the non- 
Fermi liquid behavior at the 'hot' spots is masked by the 
'cold' sections eventually restoring Fermi liquid behav- 
ior: p ~ T^. Therefore, within the MHM approach and 
in the the quantum critical region, the specific heat co- 
efficient displays divergent behavior as T — > following 
Eq. (10). The resistivity could show non- Fermi liquid 



behavior under small disorder which has been found to 
strongly influence antiferromagnetic QCP's^^. Averag- 
ing the scattering rate over the Fermi surface reduces 
the effectiveness of the Hlubina-Rice mechanism and the 
scattering from the 'hot' regions becomes effective lead- 
ing again to non-Fermi liquid behavioil^. 

The behavior of the specific heat coefficient that we 
have found around the critical point (see Fig. 10) does 



show an enhancement on approaching the QCP in con- 
sistent agreement with Eq. (10). However, we cannot ac- 



curately determine the logarithmic dependence from our 
numerical data due to the small cluster sizes reached. 



V. CONCLUSIONS AND OUTLOOK 

We have analyzed in detail the effect of geometri- 
cal frustration on charge ordering transitions realized in 
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the extended Hubbard model on the anisotropic trian- 
gular lattice, which appropriately describes the family 
of quarter-filled layered organic crystals: 0-(ET)2X. The 
model contains both onsite, U, and inter-site Coulomb re- 
pulsion terms Vp and Vc, that are taken to be isotropic, 
Vc = Vp = V. The degree of geometrical frustration 
in the electron motion is tuned through the tc/tp ratio, 
which is an important parameter controlling the experi- 
mental phase diagram. 

The zero temperature phase diagram of this model con- 
tains a homogeneous metal (HM), a pinball liquid (PL) 
and a three-fold charge ordered (SCO) phase. While the 
SCO phase occurs at sufficiently strong inter-site inter- 
actions V for any fixed U, the PL only occurs above a 
certain threshold U value, as its existence is inherently 
tied up to the strong coupling regime. On the other hand, 
the PL is found to be stabilized by increasing the geo- 
metrical frustration of the lattice. Our results do show 
that in the range of values of U/tp and V/tp appropri- 
ate to the 6'-(ET)2X materials, increasing the geometrical 
frustration of the lattice can effectively tune the system 
from a homogeneous metal (HM) with strong charge or- 
der correlations to a "pinball" liquid (PL) phase. 

The phase transitions between charge ordered and 
disordered metallic phases can display quantum criti- 
cal phenomena in close analogy with the heavy fermion 
systemJiSHIIl with the critical charge rather than the spin 
fluctuations driving the CO transition. Such type of fluc- 
tuations may be at the origin of both the anomalous 
properties in the metallic state and Cooper-pair forma- 
tion. Indeed, non-Fermi liquid behavior as well as non- 
BCS superconductivity have both been predicted and ob- 
served in quarter-flUed organic materials of the a, /3" and 
d-(ET)2^ typ6^° Such heavy fermion behavior arising 
from molecular tt electrons instead of the rf or / electrons, 
as occurs in the rare earths, may indeed find a natural 
explanation based on the properties of matter expected 
near a QCP. 

In order to establish whether quantum critical behavior 
occurs or not in the quarter-filled layered materials close 
to CO several issues could be experimentally and theoret- 
ically addressed: (i) Is there evidence for the divergence 
of the specific heat coefficient and the quasiparticle effec- 
tive mass, m* /m — !■ oo and for the collapse of the Fermi 
temperature, T* 0, near the QCP? Measurements of 
the quadratic coefficient of the resistivity approaching the 
QCP from the Fermi liquid side of the critical point can 
be useful to test the effective mass enhancement. Such 
type of experiments have been systematically performed 
in K-(DHDA-TTP)2SbF6 and (MeDH-TTP)2AsF6, by 
tuning the sy stem across the CO transition via applied 
pressure^IEIl, yielding phase diagrams similar to those of 
Fig. [9) (n) Is there non-Fermi liquid behavior of trans- 
port and thermodynamic properties in the quantum crit- 
ical regime above the QCP? What is the temperature 
dependence of the resistivity in these systems? Are there 
clear deviations from Fermi liquid behavior of the form in 



Eqs. (10) and (12)? (iii) If quantum criticality and scal- 



ing are observed in transport and thermodynamic quan- 
tities, how much of this behavior is consistent with the 
MHM predictions? Could there be a new universality 
class around the QCP in quasi-two-dimensional organic 
materials, related to the emergence of the pinball phase? 

(iv) Measurements of the Hall coefficient can be useful 
to disentangle whether the QCP is of the MHM type or 
different. In standard MHM theories, the Fermi surface 
would fold due to Bragg refiection off the density wave 
with no discontinuity in the Hall c onstant when the sys- 
tem is tuned across the QCPlSES]. 

However, as in heavy 
fermions a local type of QCP could arise in which the sys- 
tem jumps discontinously from a large Fermi surface to a 
small Fermi surface through the QCP leading to a discon- 
tinuous jump of the Hall coefficient. Since the transition 
from the HM to the PL involves localization of the 'pin' 
electrons, a transition from the large Fermi surface of 
the HM involving all carriers to a small Fermi surface in- 
volving only 'ball' itinerant electrons could indeed occur. 
Understanding how this transition takes place and the 
type of QCP observed could be resolved by Hall constant 
measurements in analogy to the heavy fermion systems. 

(v) Here we have mainly discussed transitions between 
disordered and ordered metallic phases for which MHM 
theory is meant for. An important issue to address is 
how quantum criticality is modified in transitions from 
HM to CO insulating phases? This issue can be addressed 
within the EHM studied in the present work, by allowing 
for anisotropic Coulomb interactions Vp Vc- 

As observed in Ref.l^ the superconductivity in 9- 
(ET)2X compounds, as in other polytypes, frequently 
appears near to an insulator. In such cases, the cause 
of superconductivity (SC) may not be the simple weak 
coupling BCS mechanism by the electron-phonon inter- 
action, but rather due to electronic correlations. Several 
theoretical works in the weak coupling limit have been 
performed in order to examine the possible mechanism 
for the onset of SC in proximity to the CO phaseP^^. 
Unconventional SC of the /-wave type has been en- 
countered on the anisotropic triangular lattice with the 
model parameters: tc = and Vc — Vp t^^ with 
U = lOtp and mediated by the charge fiuctuations. This 
/-wave pairing symmetry is the analogous to the d^y- 
wave pairing found in proximity to the checkerboard CO 
on the square lattice^'*. It would be interesting to search 
for unconventional SC around the QCP found for other 
tc/tp ^ ratios and U values both in the SCO and PL 
type QCP. Based on the results of the present work, it 
can be conjectured that the anomalous properties and 
unconventional superconductivity observed in X=l3 com- 
pounds maybe related to the proximity to the strong cou- 
pling PL phase since the onsite Coulomb repulsion energy 
is significant: U = (15 — 20)tp. This could be tested by 
applying uniaxial pressure on the 0-(ET)2l3 crystals. 

Other materials such as the rare-earth nickelates, 
AgNi02 , do show three- fold CO transitions similar to the 
one discussed her^SH although the origin of the CO tran- 
sition may be non-Coulomb in origin since AgNi02 has 
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a complex multiorbital structure^*' and other effects such 
as Hunds couphng and crystal fields can play a relevant 
role. Thus, the organic materials of the 6'-(ET)2X type 
appear to be ideal candidates to single out the effects on 
electronic properties of metals close to a charge-driven 
quantum phase transition mediated solely by the offsite 
Coulomb repulsion between electrons at different sites. 



Appendix A: RPA results at weak U 

In the random phase approximation (RPA) , the insta- 
bility of the homogeneous metal is signalled by a diver- 
gence of the charge susceptibilitjff^ 



x(q) 



xo(q) 



l + [;7/2 + y(q)]xo(q) 



(Al) 



at a given wavevector. Here Xo(q) is the non-interacting 
susceptibility of the lattice, ^(q) is the interaction po- 
tential in Fourier space and U is the onsite repulsion. 
For isotropic n.n interactions on the triangular lattice we 
have 



Viq) = 2V{cosiq,) + 

+ cos[(g^ + V3qy)/2] 



(A2) 



cos 



iq., - V3qy)/2]}. 



An instability occurs when the denominator in Eq. ( Al ) 
vanishes, which requires — V^(q) = Xo ^(q) + U/2. In 
principle the above equation can describe both charge 
ordering, driven by the Coulomb iteraction — V^(q) that 
is maximum at the six equivalent threefold wavevectors 
Q = (±27r/3,±27r/y3), (±47r/3,0) (blue dot in Fig. 
[TT^), and a charge density wave (CDW) induced by a 
large xo • The evolution of the free-electron susceptibility 
with frustration is illustrated in Fig. |11| 3. 

A CDW instability occurs for tc/tp = —1 due to a 
perfect nesting between parallel segments of the Fermi 
surface (red lines in Fig. [TT^), at wavevector = 
(tt, 7r/-\/3) (red dot, corresponding to the M-point), but 
is washed out at U ftp > 7, where the threefold order is 
always favored. An instability also appears to compete 
with the 3C0 in the case tc/tp = 1 for [/ < 2tp, at a 
wavevector Qi — (2.62,2.72) (white dot). Such vector 
lies at the intersect between a circle of radius |Qi| — 2k p 
and the boundary of the Brillouin zone. It represents 
a compromise between a Fermi surface instability and a 
genuine charge ordering, as it benefits from both a large 
Xo and a large —V{q). For values of the frustration ratio 
\tc/tp\ < 0.9, the RPA predicts that the 3C0 transition 
is dominant for all C/ > 0. 

The critical coupling V^co at the threefold instability 
is shown as a function of the geometrical frustration in 
Fig. [TT| 3. The RPA predicts an increase of V^co with 
tc/tp which originates from a decrease of the density of 
states at the Fermi leveP^. We see that even at a rel- 
atively low value oi U = 5tp, the RPA result does not 
agree with the exact diagonalization data, showing an 
opposite trend for positive tc/tp ratios. In the ED, the tc 




K' M' 



M K 



K' M 



FIG. 11: (Color online) (a) Brillouin zone of the triangular 
lattice (black), reciprocal lattice vectors Gi — (27r, — 27r/\/3) 
and G2 = (0,47r/\/3) (arrows), and Fermi surface: from light 
blue to dark blue, tc/tp = 1,0.5,0,-0.5; red, tc/tp = —1 
(same color code as in Fig. [2|. The blue and red dots are 
respectively the threefold wavevector Q = K = (Ivr/S, 0) and 
the nesting wavevector Qp — M. — {n, 7r/\/3) connecting flat 
segments of the Fermi surface at tc/tp = — 1 (red), while the 
white dot is the wavevector associated with the predominant 
instability for tc/tp = 1; (b) the critical value of the intersite 
Coulomb interaction Vsco for U = 5ip, comparing the RPA 
and ED results, (c) Free-electron susceptibility xo along sym- 
metry lines of the Brillouin zone, for different values of tc/tp. 



dependence is governed by the stabilization of the pinball 
liquid phase, that is not captured by the weak coupling 
RPA argument. 



Appendix B: Mean-field potential energy in the CO 

state 



The Hartree expression for the potential energy per 
site in a charge ordered state with three-fold symmetry 
reads: 



E 



U 



H 



(nA-triAi + riB^nBi + nc^nci) 



+ V{nAnB + riAUc + riBnc). 



(Bl) 



We take = tiai = ub^ = ubi = 1, nc^ = 1/2, nci = 
corresponding to the Hartree solution for the 3C0 state 



and nA-[ = riBf = 3/4, uai = 
for the pinball liquid phase 



1, 



net 



1, 







^see Figl4|. Inserting these 
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values in the preceding expression we find: 



E' 



(M) 



E 



H 

(PL) 
H 
(3CO) 
H 



E 



(9/16)[/ + (27/4)y, HM (B2) 
(l/2)C/+(13/2 + l/16)y, PL (B3) 
(2/3)C/ + 6y, 3C0 (B4) 



The potential energy calculated by ED closely follows 
the V dependence predicted by the above mean-field 

equations (not shown). The transition from the pinball 
state to the three-fold CO state is well captured by the 



mean- field analysis: the ED data closely follow the value 
Vc = U/3 obtained by equating E^^^^ = E^^'^^^ . 
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degree of frustration in the electron motion via the parame- 
ter tc/tp. Previous studies of the EHM on the square lattice 
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dependent of tc/tp in the small A^c = 12 cluster used which 
coincides with the tc = kinetic energy of the extended 
tight-binding model of the lattice. At > 5t dcgenoracios 
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